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We present a theoretical treatment for the dissipative two-resonator circuit quantum electrody- 
namics setup referred to as quantum switch. There, switchable coupling between two superconduct- 
ing resonators is mediated by a superconducting qubit operating in the dispersive regime, where 
the qubit transition frequency is far detuned from those of the resonators. We derive an effective 
Hamiltonian for the quantum switch beyond the rotating wave approximation and provide a detailed 
study of the dissipative dynamics. As a central finding, we derive analytically how the qubit affects 
the quantum switch even if the qubit has no dynamics, and we estimate the strength of this influ- 
ence. The analytical results are corroborated by numerical calculations, where coherent oscillations 
between the resonators, the decay of coherent and Fock states, and the decay of resonator-resonator 
entanglement are studied. Finally, we suggest an experimental protocol for extracting the damping 
constants of qubit and resonators by measuring the quadratures of the resonator fields. 

PACS numbers: 84.30. Bv, 03.67.Lx, 32.60.-|-i, 42.50.Pq 
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INTRODUCTION 



Circuit quantum clcctrodynamicsl"'^ (QED) is the 
solid-state analog of quantum-optical cavity QED.— 
While in the latter natural atoms are coupled to 3D- 
cavities, the former is based on superconducting quantum 
circuits and the roles of the atoms and the cavities are 
played by qubili^^— and microwave resonator circuitai^rJ^, 
respectively. In fundamental research, circuit QED ar- 
chitectures have proved to be valuable for implementing 
quantum optics on a chip, for which a rich toolbox has 
been developedJ^"— These experiments were based on 
a single qubit coupled to a single resonator. With ap- 
plications for quantum information processing in mind, 
an extension to multiple qubits seems naturali^"— More 
recently, the potential of using multiple resonators has 
been pointed out by several authorsi^i"— 

Under opportune circumstances, in a two- resonator 
circuit QED setup, a superconducting qubit acts as a 
quantum switch between two superconducting on-chip 
rcsonatorsj^ To this end, the qubit must be detuned 
from both resonators. The resulting effective Hamilto- 
nian describes a resonator-resonator interaction whose 
coefficient has two contributions. The first contribution 
depends on the qubit state and the qubit-resonator de- 
tuning and can have a positive or negative sign. The 
second one has a definite sign and stems from the fact 
that qubit and resonators arc not point-like objects but 
extended circuits. Provided that the qubit always is 
in a suitable energy eigenstate, the switch is turned off 



when both terms are balanced and turned on otherwise. 
Beyond this simple protocol, the "quantumness" of the 
setup can be exploited by bringing the qubit into a super- 
position state with the resonators. This allows for gener- 
ating bi- and tripartite entanglement or Schrodinger cat 
states. 

In a real experiment, one expects the operation of the 
two-resonator circuit QED setup to be affected by the 
various decoherence rates of qubit and resonators. Since 
most implementations of superconducting qubits can be 
tuned by external parameters, those rates depend not 
only on the qubit type, but also on the operating point. 
So far^"— they have been in the range of approximately 
1-200 MHz. In the case of the quantum switch, only qual- 
itative estimates on the effect of qubit dephasing exist 
However, a detailed quantitative understanding of the 
possible effects stemming from the various existing de- 
coherence channels is indispensable for successful exper- 
imental implementation. In particular, it is essential to 
analyze the effect of qubit decoherence sources on the 
coupled resonator pair. Hence, in this work, we develop 
a complete dissipative theory for a circuit QED setup 
consisting of two resonators both dispersively coupled to 
a single qubit. As a central result, we demonstrate that 
qubit relaxation affects the resonators in second disper- 
sive order, whereas dephasing becomes an issue only in 
fourth dispersive order. 

The paper is structured as follows. In Sec. [TTl we in- 
troduce the system Hamiltonian and add the baths caus- 
ing dissipation to the system. We model the bath influ- 
ence using a Bloch-Redficld quantum master equation. 
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Next, in Sec. IIIIBl wc derive an effective system Hamil- 
tonian beyond the rotating-wave approximation (RWA). 
We show that this extension results in quantitative, but 
not quahtative changes compared to a treatment within 
RWA. Furthermore, in Sec. IIII C[ we derive a simplified 
effective quantum master equation suitable for analytical 
treatments. In particular, we use the latter result to com- 
pute an explicit expression for the influence of the qubit 
dissipation channel on the two-resonator system. Sec. lIVI 
contains numerical results for various prototypical oper- 
ation modes of the quantum switch setup. These include 
coherent oscillations between the resonators, the decay 
of Fock and coherent states, and the decay of resonator- 
resonator entanglement. We show that the agreement 
with the analytical results obtained by means of the ef- 
fective quantum master equation of Sec. IIII Cl is excellent. 
Most importantly, we show that qubit dissipation affects 
the switch only in second dispersive order. Finally, in 
Sec. El we suggest a protocol to extract the damping 
constants of the system by measuring the field modes of 
the resonators. The appendices contain technical details 
about the calculations presented in this article. 

II. DISSIPATIVE TWO-RESONATOR CIRCUIT 

QED 

We introduce a dissipative description for a circuit 
QED architecture consisting of two on-chip microwave 
resonators that are simultaneously coupled to one su- 
perconducting qubit. This setup is sketched in Fig. [T] 
We emphasize that our formalism is general in the sense 
that qubit and resonators can be based on any suitable 
quantum circuits. However, whenever we need to give 
numbers, we assume a persistent-current flux qubilj^i^ 
coupled to two transmission line resonators henceforth. 

A. System-bath model 

First, we write down the two-resonator circuit QED 
Hamiltonian. A detailed derivation is given in Ref. H^- 
The natural reference frame is the laboratory basis, the 
physical basis of circuits and fields. 

n' =—0-' + ^cf' + hVlAa'^a + hVtBb'^b 
2 2 

+ ftG(a + at)(5 + &t) (1) 

+ hgA(T',{a + a^) + hgncr'^ib + b^) . 

The first two terms in the first line of the above Hamil- 
tonian represent the qubit in terms of the standard Pauli 
operators a'^ and cr^. The controllable energy bias is he, 
and HSq denotes the minimum level splitting. In the par- 
ticular case of a flux qubit j^i^ Mg is the tunnel split- 
ting, and the energy bias He = 2iq($^^— $o/2) can be 
tuned by an externally applied flux The quanti- 

ties iq and $o — h/2e denote the qubit persistent cur- 
rent and the magnetic flux quantum, respectively. When 
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FIG. 1: (Color online) Sketch of the two-resonator circuit 
QED system under analysis, including schematically the in- 
teraction with bosonic heat baths (blue boxes with oscillator 
potentials). Both microstrip or coplanar waveguides could be 
employed as resonators (red lines). As in the text, the cou- 
pling qubit is exemplarily depicted as persistent-current flux 
qubit (green loop). The system is coupled to external circuits 
via coupling capacitors. The decay rates 7 and 7^ are 

defined in Sec. HTCl 

"^x"" ~ 'I'o/S or, equivalently, £ = 0, the qubit is said 
to be biased at its degeneracy or optimal point, where it 
is protected from low-frequency noise to first order. The 
last two terms in the first line of Eq. ([T]) represent the two 
resonators with frequencies and fls- Here, a,b and 
a' , b^ are the annihilation and creation operators of the 
modes in resonators A and B, respectively. The second 
line of Eq. ([T]) describes the geometric coupling between 
the resonators, which is due to the fact that we are deal- 
ing with circuits. The coupling coefficient G contains 
contributions both from a direct coupling and an inter- 
action that is mediated by the qubit circuit. Finally, the 
third line of Eq. ([T]) describes the qubit-rcsonator cou- 
pling terms with coefficients gA and gs- As explained 
in Sec. IIII B| they give rise to a "dynamical" resonator- 
resonator coupling under appropriate conditions. 

In a real experimental scenario, the two-resonator cir- 
cuit is unavoidably coupled to an external circuit that 
is characterized by an impedance Z{uj). In a quantum 
mechanical description, this impedance can be modeled 
by coupling the circuit bilinearly to the modes of an elec- 
tromagnetic environment consisting of an infinite set of 
harmonic oscillators Following this route, we obtain 
a Caldeira-Leggett-type system-bath Hamiltonian,—"— 

+ EEk(4a. + ^) (2) 

The indices ji G {A,B,x,z} label the system-bath cou- 
pling operators with respect to the different reservoirs 
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the system is coupled to. In detail, 

Qa = (a + a^) , Q: 
QB = {b + b^), Q. 



a'. 



(3) 



The coupling coefficients represent the interaction be- 
tween the system and the different bath modes with fre- 
quencies cj^, which are described by the bosonic annihila- 
tion and creation operators dj^^,dj ^. Within the scope 
of this paper, we consider the noise sources to be un- 
correlated. This is justified since the different types of 
noise are caused by fluctuations of distinct nature. In 
other words, we assume that the baths are independent, 
= SijS^„. We find it noteworthy to mention 
that for n S {a;, z} the coefficients Cj depend on the spe- 
cific implementation of the qubit . For a flux qubit, the 
dominant noise source is believed to be flux noise,—"— 
which couples to the circuit via the z-axis in the labora- 
tory frame. 

In order to get more physical insight, we rotate H' into 
the qubit energy eigenbasis^ {|g),|e)}, where \g) and 
|e) denote the flux-dependent qubit ground and excited 
state, respectively. Using the redefined Pauli operators 



smf (T, 



f^x = |5>(e| + |e)(.g| ^cos9( 
<^z = |e) (e| - \g) {g\ = sin 9(j'^+ cos 9 a'^ 



(4) 



we obtain 



n 



'qb. 



+ hG{a + a^){b + b'^) 

-\- hgA (cos 9 Gz — sin 9 a,^){a + a)) 

+ hgB{cos 9 Oz— sin 9 (Tx){b + 6^) . 



(5) 



The flux-dependence is now encoded in the qubit energy 
level splitting fiwqb = ^(^q -I-£^)^/^ and the mixing angle 
9 = arctan(5Q/£). The qubit-batli coupling operators are 
rewritten as 



Qx 



(^x 



- COS V Ux 
COS 9 Oz 



- sm a Oz 
sin 9 (Tr 



(6) 



They are deflned along the rotated axes determined by 
the tunneling matrix element Ti^q in a^-direction, and 
the energy bias fe in cr^-direction. The system-bath in- 
teraction is fully characterized by the spectral densities 



(7) 



In the case where decoherence is mainly caused by exter- 
nal circuitry, the spectral densities are proportional to 
the real part of the impedances Re[Z^(w)]. In general, 
internal loss mechanisms are also relevant in supercon- 
ducting resonators at low powers and low temperatures. 
They often originate from fluctuators on the resonator 
surface, which are usually modeled as two-level systems. 



Thus, we interpret the Jfj,{uj) in an effective sense in 
that they include both the effects of external circuitry 
and internal losses. Our effective description does not 
cover the so-called excess phase noise though, i. e. low- 
frequency fluctuations in the resonator frequency itself, 
which originate from the surface fiuctuators as well. As 
it was pointed out and investigated experimentally ; ^'^i^^ 
this leads to resonator dephasing. While such effects are 
not included in our modeling of decoherence, we cannot 
ensure that they will only be of minor importance with 
respect to operating the two-resonator setup (see below 
in Sec. IIII A)) . In most experimental situations, however, 
decoherence is predominantly governed by external res- 
onator losses. The corresponding external quality factor 
is characterized by the coupling capacitors to external 
circuitry. We note that resonator dephasing was not re- 
ported to play a major role in recent circuit QED experi- 
ments done with comparable resonators. In any case, the 
role of non- vanishing excess phase noise requires a sepa- 
rate, more detailed treatment with respect to an intended 
experimental realization of our setup. 



B. Bloch-Redfield quantum master equation 

The dissipative dynamics of the qubit-two-resonator 
system is obtained by tracing out the bath degrees of 
freedom of the total density operator gitot associated with 
the transformed system-bath Hamiltonian, 



(8) 



where the qubit-bath coupling operators Qx and Qz are 
now written in the qubit eigenbasis according to Eq. ([6|) . 
For weak system-bath interaction, the baths can be elim- 
inated within Bloch-Redfield theory^i^ as follows: As- 
suming that the baths are initially in thermal equilib- 
rium at temperatures and not correlated with the 
system state g, the total system-bath state can be writ- 
ten as gtotocg (g) exp(- V ^']d],^d.^^/kY>,T^). Then, 
one can derive within perturbation theory the quantum 
master equation for the reduced system density operator 
Q = Trbath[etot]- This procedure yields 



e(i) = --[H,f?(i)] 



(9) 



1 1"^ 



h.c. 



The environment correlation functions K^{t) are given 

by 
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coth 



fvjj 



cos LOT— I smujT 



(10) 
The 



where, J^(w) are the spectral densities dZ)) 
Heisenberg operators Qiiir) = uI{t)Q f^UQ{T) are 
constructed via the system propagator Uo{t) = 
T[exp{— /J" diH(t)}]. Here, the time ordering oper- 
ator T is only required for an explicitly timc-depcndcnt 
system Hamiltonian. 

We note that Eq. ([9)) is based on a Born-Markov ap- 
proximation, since the bath correlation functions are sup- 
posed to decay sufficiently fast as compared to typical 
timescales of intrinsic system evolution. Thus, it was 
appropriate to extend the integral in Eq. ([9]) to infin- 
ity. Consistently, we assume Ohmic spectral densities in 
the correlation functions of Eq. pil)) . modeling Z{ijj) as 
an effective resistance. However, this restriction is only 
necessary in the low-frequency region of the qubit envi- 
ronments. There, we assume 

J^_,{uj) ^ a^u /.t e w<a;qb (11) 

and the coefficients represent the dimensionless damp- 
ing strengths. As we will see later, in the high-frequency 
regime, we are interested only in infinitely small in- 
tervals around frequencies such as Wqb, ^a, and ^Ib- 
Hence, the Born-Markov approximation remains justi- 
fied by expanding J^{lo) to first order in these intervals. 
In this way, the only remaining restriction is that J^{uj) 
is a smooth function around the frequencies of relevance. 
Within the scope of this work, we shall consider Eqs. 
pT|) as a full description of the influence of dissipation 
and decoherence on the two-resonator setup. 

This reasoning excludes in particular l//-noise, 
which affects the phase coherence of superconducting 
qubits^i^i^ due to its high impact at low frequencies. 
One typically describes l//-noisc by calculating the ac- 
cumulated random phase as a function of time for specific 
experimental protocolsi^^illiM However, as shown below, 
we expect the effect of qubit dephasing to be suppressed 
even more than relaxation effects in the setup described 
here. 

Since the quantum master equation ^ is non-trivial 
with respect to analytical treatment, we only use it for 
numerical purposes. However, in Sec. Illli we derive a 
simplified effective quantum master equation in the dis- 
persive regime, which will allow for analytic insight into 
the dissipative behavior of the two- resonator circuit. 



work in the qubit energy eigcnframc and refer to what 
is commonly called energy relaxation and pure dephas- 
ing. The former describes bath-induced level transitions, 
while the latter accounts for the pure loss of phase co- 
herence without a change of the system energy. In order 
to define the decay rates corresponding to these two pro- 
cesses, we first review the dissipation mechanisms of the 
qubit alone. To this end, we derive the quantum master 
equation describing the evolution of the reduced qubit 
density matrix Qqh for a single qubit associated with the 
Hamiltonian 'Hqh = ftaiqbcrz/S in the energy eigenbasis. 
Using the formalism detailed in App. \^ we find 

+ 7(^qb) (^cr^ gqh{t)<j'^ - ^[cr^cr",£'qb(i)]_ 

(12) 

where i5]_|_ = AB+BA denotes the anti-commutator 
between the operators A and B. The dissipator in the 
third line of Eq. (|12p does not affect the populations of 
the qubit eigenstates, but only accounts for the decay of 
the off-diagonal elements of the density operator. Thus, 
the rate 70(^^0) can be associated with pure dephas- 
ing. The dissipator in the second line of Eq. ((T^ induces 
transitions between the qubit eigenstates, hence 7(wqb) 
characterizes relaxation. Assuming an overall tempera- 
ture T = T^=T^, and following Eq. (|A7| and Eq. (jAQ]) . 
the qubit energy relaxation rate 7(a;qb) and pure dephas- 
ing rate 70(0;— >-0) are obtained as 



7(Wqb) = Jxi^qh) COS^ ^ + Jzi^qh) siu 



70(^^0) 



cos^ 9 + ax sin^ 



(13) 
(14) 



Equations p3|) and ([T4|) link the physical system-bath in- 
teractions quantified in the laboratory frame to the pure 
bit-flip and dephasing mechanisms relevant in the qubit 
eigenbasis. Moreover, they highlight the dependence of 
the pure decay rates on the applied flux in terms of the 
mixing angle 6. In particular, for a fiux qubit, flux noise 
can be responsible for both relaxation and dephasing. We 
emphasize that, in this special scenario, Ja;(aj) = and 
Jz{io) ^0, and Eq. (fT^ is consistent with results from 
other works i^MLda 



III. ANALYTICAL TREATMENT OF 
DECOHERENCE IN THE DISPERSIVE LIMIT 



C. Qubit decay rates 

So far, we have modeled the coupling of the qubit to 
the baths in the laboratory frame. In this way, we can in- 
clude the relevant noise channels for any particular qubit 
architecture into our formalism easily. However, with re- 
gard to physical understanding, it is more favorable to 



In the setup of Fig. [1] the qubit can mediate a control- 
lable coupling between the two resonators, i. e., it can act 
as a quantum switch between them. In this section, we 
review the quantum switch Hamiltonian of Ref. [l^ and 
extend it beyond the rotating-wavc approximation. Fur- 
thermore, wc derive an effective quantum master equa- 
tion which allows us to understand by purely analytical 
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arguments that the quantum switch is affected by the 
qubit dissipation only in second (relaxation) and fourth 
order (dephasing), respectively. 



A. Dispersive Hamiltonian within the rotating 
wave approximation: The quantum switch 

In order to function as a quantum switch, the two- 
resonator circuit must be operated in the dispersive limit, 
where the qubit-resonator detuning A is large as com- 
pared to the qubit-oscillator coupling, 

g < A , A^LUqb-n, (15) 

and the parameter Aa is necessarily small: 



mg 



Aa 



gsinO 



|Aa| « 1 



(16) 



Here and henceforth, we confine ourselves to symmetric 
setups with n — ^{a.b} ^-nd g = g{A.B}- This is not 
expected to be a serious restriction in practice, though.— 
In the dispersive limit determined by Eq. (|15|). the 
Hamiltonian of Eq. ([Sj can be diagonalized approxi- 
mately. To this end, it is first simplified with a rotating 
wave approximation as follows. Writing = + cr^ 
with the fermionic raising and lowering operators (t+ = 
|e)(g| and = \g){e\, one can move to the interac- 
tion picture with respect to the uncoupled Hamiltonian. 
Then, the coupling operators a~^a, (j^a\ a^b, and a^b^ 
oscillate with the phase factors exp[±iAi], whereas cr^a, 
a~^a\ a~b, and a'^b^ oscillate with exp[±iEt], where 



fl + UJqh ■ 



(17) 



Close to resonance, the resonator-qubit detuning is small 
and, consequently, |A| S. Thus, the former set of 
operators oscillate slowly, whereas the latter exhibit fast 
"counter-rotating" oscillations. For sufficiently weak cou- 
pling g <ti minjwqb, ^2}, one can separate timescales and 
average the counter-rotating terms to zero. In this way, 
the first-order interaction Hamiltonian between qubit 
and resonators is Jaynes-Cummings-like^ and we de- 
scribe our system with 



n 



RWA 



hilb^'b 



- hAXA{(7+a + (j-a^ +a+b + (T-b'^) (^^^ 
+ hG{a^b + ab^) . 

In a second step, we apply the transformation ^^WA _ 
exp(— Aa2^), where 



V = <T^a^ - a+a + a^b^ - a+b . 



(19) 



Finally, one truncates the transformed Hamiltonian 
■HRWA ^ ul^^n^^^UBWA to second order in Aa, yicld- 



K 



RWA 



huj, 



qb 



:nr2(a^a + 6^6 + 1) + 

+ hAX\<7^ {a^a + b'^b+ l) 
+ Sgs^^^(a6t + at6). 



(20) 



Here, the first line describes qubit and resonators, the 
second ac- Stark/ Zeeman and Lamb-shifts, and the third 
an effective coupling between the two resonators. The 
corresponding coupling operator is 



'RWA 
i/SW 



G 



(21) 



A remarkable feature of the Hamiltonian of Eq. (|20l) is 
that it commutes with CTz, i.e., [^JJff^^, fz] = 0. Conse- 
quently, the qubit state will not change during the uni- 
tary evolution of the system. When the qubit is prepared 
in a suitable eigenstate, it can be traced out. Throughout 
this work, we consider the qubit to be initially prepared 
in its ground state \g){g\- Then, gsw simplifies to the 
resonator-resonator coupling constant 



9sw^ — G — A^A . 



(22) 



In this case, the Hamiltonian of Eq. ([20|) describes two 
coupled harmonic oscillators. By means of either adia- 
batic or oscillating external fiux signals, the qubit can be 
tuned such that the interaction between the resonators 
is cither switched on (|gf^^| ^ 0) or off {\gf^^\ = 0). 
This feature is referred to as the switch-setting condition. 
With the help of specific protocols, it can be utilized to 
create entangled states out of initial bi-resonator product 
states. 

We note that the effective coupling between both res- 
onators can be interpreted as a beam-splitter interac- 
tion. A comparable quantum optical setup was proposed 
in Ref. [sol There, an atom passing through a cavity 
serves to create entanglement between two optical fields 
inside the cavity. That system is described by an effective 
Hamiltonian quite analogous to Eq. ([20)) . 

For this work the "adiabatic" shift protocol is of par- 
ticular relevance. There, parameters are initially chosen 
so as to fulfill the switch-setting condition when the qubit 
is in \g). Then, the resonator-resonator interaction can 
be turned on by adiabatically varying the flux bias. Ex- 
periments have shown that a flux change slow enough to 
avoid significant population of the excited state can be 
realized easily even in pulsed setups 

Regarding the influence of a dissipative environment 
on the quantum switch, we already note that, for suit- 
able switching protocols and at sufficiently low tempera- 
tures, the qubit energy relaxation and dephasing will not 
affect the operation of the switch in first order. As one 
of the main results of this work, we give analytic argu- 
ments to put this statement on firm theoretical footings 
in Sec. IhTCI 
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B. Dispersive Hamiltonian beyond the 
rotating-wave approximation 

In the process of deriving Eq. in the previous sec- 
tion, a rotating wave approximation is apphed to the 
Hamiltonian ([S]) at the level of first-order in the qubit- 
oscillator coupling. However, it has recently been re- 
vealed that neglecting the counter-rotating terms may 
lead to inaccuracies.— Especially in the case of far de- 
tuning described by Eq. (|15p . the rotating-wave approxi- 
mation causes noticeable deviations from results obtained 
numerically from the full Hamiltonian of Eq. ^ for typ- 
ical parameters. Nevertheless, the effective, dispersive 
Hamiltonian can be obtained by means of the unitary 
transformation 

= exp ( - XaV - Ae5 - XnW) . (23) 

Here, 

S = a-a- +a-b- a+b"^ (24) 

W = a,ia~a^) + a,ib-b^), (25) 

and the corresponding coefficients are 

A.^^, |A.|«1 (26) 

Ao^^, |Ao|«l. (27) 

The above inequalities allow us to discard terms of or- 
ders higher than A^^ ^ when computing the effective 

second-order Hamiltonian HeS = W^WA. In this case, 

U,^ = hn{a^a + b'^b+l) + ^a, 

+ ft(AAi + SA|) (7^0^0 + 6^6 + 1) (28) 
+ figsw{ab^ + a^b) 

becomes diagonal. In the above equation, we use the 
qubit-state-dependent resonator-resonator coupling op- 
erator 




0.495 0.496 0.497 0.498 

FIG. 2: (Color online) Switch setting coefficient gsw (qubit 
in the ground state) for the RWA (blue dashed line) and non- 
RWA case, Eq. (|30|) (red dashed line), compared to numerical 
data obtained from the diagonalization of Hamiltonian ^ 
(black solid line). All quantities are plotted against the qubit 
flux bias <i>J!"^. Parameters are Q/2-k = 3.5 GHz, 5q/2tt = 
4 GHz, G/2tt = 0.0022 GHz and g/2TT = 0.24 GHz. 



However, in HoS, terms of the order X^^^^^^, such 

as, e.g., cr+(a(^^')^ or a^a^a, need to be canceled with 
a rotating-wave argument. We emphasize that this 
rotating-wave approximation in second-order in Heff still 
allows for an accurate description of our system in the 
dispersive regime, whereas a rotating-wave approxima- 
tion in the first-order Hamiltonian H docs not. Fol- 
lowing the same reasoning, we may also neglect terms 
(X G/A, G/S,... ^ 1. The effective Hamiltonian Hcs 
of Eq. (pS)) has the same structure as its rotating-wave 
counterpart H^^^ of Eq. (PU)) . However, there is one 
important quantitative difference: the detuning depen- 
dence A|^A of the coupling coefficients is replaced by 
the expression {Xj^A + X^J^). In particular, the effec- 
tive resonator-resonator coupling constant is 



gsw^G- {XlA + A|E) (30) 



5sw = G+(AiA + A|S)a,. (29) 

At this point, we can gain insight about the effect of 
the transformation U [Eq. (|23p ] on the Hamiltonian. We 
first note that, when applying the rotating-wave approx- 
imation to "H, only the exponent V is required to pro- 
duce a diagonal second-order Hamiltonian, cf. Scc. lIII Al 
Beyond this simple scenario;^ the exponent S cancels 
the first-order counter-rotating terms of Ji. Further- 
more, the polaron transformation represented by the ex- 
ponent W must be applied to eliminate off-diagonal in- 
teraction terms such as g cos 9a ^ {a -I- a^), which cause 
qubit-state-dependent energy shifts of the oscillator co- 
ordinates when the qubit is biased away from its degen- 
eracy point. 



for the qubit being in its ground state. The effect of the 
counter-rotating terms is visualized in Fig. O There, we 
compare ffg^^ s-i^d 5sw to the numerically exact cou- 
pling coefficient for adequate parameters. Obviously, in 
contrast to ffgw^j the agreement is excellent for gsw- 
This finding once more confirms the necessity to in- 
clude counter-rotating terms of first order in the qubit- 
oscillator coupling in the full system Hamiltonian. It 
also confirms the validity of the rotating-wave approxi- 
mation in second order of A{A.s,a} applied to the dis- 
persive Hamiltonian. We also illustrate the importance 
of the non-RWA features below, where we develop a dis- 
sipativc description of the quantum switch Hamiltonian 
coupled to different reservoirs. 
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C. Effective master equation for the quantum 
switch setup 

In this section, we analytically investigate the dissipa- 
tive behavior of the two-resonator-qubit system. To this 
end, we derive an effective quantum master equation for 
the reduced density matrix of our system in the dispersive 
limiti^i^ In particular, we study the additional dissipa- 
tion imposed on the resonators due to the presence of the 
qubit. 

In principle, we combine the procedure explained in 
Sec. IIIIBI with that of Sec. Ill CI and apply it to the 
system-bath Hamiltonian Titot of Eq. ([5]) , which includes 
all counter-rotating terms. First, we compute the to- 
tal dispersive Hamiltonian 'Htot,Gff = U^Htod^ using the 
transformation (j23p and truncate it to second order with 
respect to the parameters X{A.^,n}i as described above. 
During this procedure, we obtain the effective system- 
bath coupling operators 



(31) 



The explicit expressions for these effective coupling op- 
erators are given in App. [B] In the next step, we derive 
the effective quantum master equation following the lines 
of Refs. I54li55l . While the interested reader can find the 
details in App. [C] we give a short summary of the most 
important steps in the following. Motivated by the usual 
experimental conditions in circuit QED, we assume an 
equal temperature for all baths, T = T'{a;,z,y!i,B}j and con- 
fine ourselves to the low-temperature regime k^T /h <C 
minjwqb. A, f2, E}. Consequently, we neglect all contri- 
butions to the dissipative system dynamics that describe 
energy absorption from the baths. 

Using Eq. Q as a starting point, we move first to 
an interaction picture with respect to the uncoupled 
qubit and resonators and insert the spectral decompo- 
sitions of the effective coupling operators. In the follow- 
ing, we perform a semi-secular approximation)^ To this 
end, we dismiss terms that evolve rapidly compared to 
the time evolution of the system state, i. e. on system 
timescales {ri,ajqb, A, E}"-"^. On the contrary, we keep 
those terms that oscillate slowly at frequencies such as 
A^A, A|;E, AqQ. We emphasize that our result goes be- 
yond the standard Lindblad master equation, where one 
would perform a full secular approximation, dismissing 
all oscillating contributions. In this way, we obtain the ef- 
fective quantum master equation for the reduced system 
state, Eq. (IC17I) . There, we assume a^^ ^,^ <C hoj^h/kBT 
in order not to violate the Markov approximation. 

In order to gain physical insight into the influence of 
dissipation on the quantum switch setup, we can simplify 
the complicated effective master equation of Eq. (|C17p . 
In the dispersive regime, the qubit mediates part of the 
coupling between the resonators by exchanging virtual, 
but not real excitations with them. In particular, as dis- 
cussed in Sec. IIIIBI the switch can be operated in a way 
that the qubit is initially prepared in the ground state 
and remains there during the whole time evolution, as it 



cannot suffer from further decoherence. In this scenario, 
the qubit degrees of freedom can be traced out and the 
reduced Hamiltonian of the coupled resonators becomes 

= hn{a^a + b^b+l) 

-ft(AiA + A|E) (a^a + 6^&+l) (32) 
+ %sw(a&^ + a^b) . 

With the dissipator acting on an operator X in the 

product Hilbert space of the resonators, 

V[X]g,,, = Xg,,,X^ - ^ [X^X, ecav]+ , (33) 

we can write down the effective Lindblad-typc quantum 
master equation for the reduced state pcav of the two 
coupled oscillators up to second order in Aa and As, 



o/cff 

'*'cav 7 t/cav 



-I- KaT) [a] Pcav + KbT) [b] Qc 
+ KqbP [a -1-6] Pcav • 



(34) 



The above equation reveals the relevant processes gov- 
erning the dissipative behavior of the quantum switch. 
The dissipators 2? [a] and 2?[fe] represent the independent 
decay channels due to the individual environments of 
the resonators A and B, respectively. The correspond- 
ing decay rates are the inverse lifetimes of the uncou- 
pled resonators, ka = JAi^) and kb = Jsi^)- These 
rates may incorporate the combined effects of internal 
and external loss mechanisms, according to the discus- 
sion in Sec. Ill Al In addition to these contributions, the 
qubit introduces extra dissipation on the resonators via 
the dissipator I?[a -I- b]. The appearance of the center- 
of-mass coordinate a+b of the two-resonator system in 
the dissipator originates from the system Hamiltonian of 
Eq. ([5]) , where the qubit couples to the resonator "center 
of mass" coordinate, i. e. the interaction is proportional 
to ax{a + b + a'^ + b'^) and az{a + b + a^ + b'^), respectively. 
The qubit-induced damping rate is 



Kqb = (Aa+As)^ {J.^{n) cos^ + Jz{n) sir 



(35) 



where 7(ri) is the rate defined in Eq. (|A7[) for the bare 
qubit. In the expressions for ka, kb, and Kqb, the spec- 
tral densities J{a,b,x,z}{^) are required to be smooth 
functions at a; = in order that Ohmic behavior can 
be assumed locally. The reasoning is the same as the one 
presented in App. [Xj 

The qubit-induced damping rate of the two-resonator 
system, Kqb of Eq. ([55]). constitutes one central result 
of this work and has several remarkable features. First 
of all, we note that 7(ri) has the same functional de- 
pendence on the qubit mixing angle 6 as the relax- 
ation rate 7(a;qb) of the bare qubit, Eq. (|13p . However, 
J{A,B,x,z}i^) and the corresponding J{A,B,x,z}ii^qh) are 
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not necessarily equal, thus the values of both rates are 
different in general. Second, the rate Kqh is of second or- 
der in Aa and As because the qubit-mediated interaction 
responsible for the effective noise channel in Eq. p4p is 
a second-order effect. This also explains the, at a first 
glance, surprising fact that the qubit induces a decay of 
the two-rcsonator system even though its excited state 
is never populated. We can understand this by recall- 
ing that the resonator-resonator interaction is mediated 
not by real, but by virtual qubit excitations, which are 
known to give rise to second-order effects. Equivalently, 
we may apply a more classical picture, which is based on 
the fact that the resonator-resonator coupling coefhcient 
f;sw of Eq. ([30|) depends on the qubit-resonator detun- 
ing. Hence, the qubit baths, which cause first-order fluc- 
tuations to the qubit level splitting, induce second-order 
fluctuations of .gsw- The latter are described by the last 
term of the master equation ([M]) . 

Remarkably, the associated decay rate Kqb is related to 
the qubit relaxation 7, whereas dephasing 7^ would en- 
ter the effective master equation, Eq. (|34p , only in fourth 
order in A{a,e} [cf. also Eq. (|Cf 7p ]. Mathematically, this 
can be understood from the structure of the dispersive 
operator (Tz^eS of Eq. (|B5|) . which couples the system to 
dephasing noise. To the order A{a,s}i this operator con- 
tains products of operators which change the populations 
of the qubit and resonators simultaneously. On the one 
hand, the term cr+a, for example, describes the excita- 
tion of the qubit together with the emission of a resonator 
photon, a process which is energetically forbidden at low 
temperatures for A = — > 0. On the other hand, 
terms such as a^a^ and a^b^ have no effect when the 
qubit remains in the ground state. By contrast, the op- 
erator ax,eS of Eq. ()B4| . which is responsible for the qubit 
energy relaxation, contains terms such as ct^ (a -I- b) of the 
order A{a.s}- These describe a resonator decay with- 
out exciting the qubit, which is energetically favorable at 
low temperatures. For this reason, the only remaining 
contribution to qubit-enhanccd decay up to second or- 
der in A{a.s} in the effective quantum master equation 
for the two resonators, Eq. ()34|) . stems from qubit relax- 
ation. The fourth-order contribution to the dephasing is 
related to the appearance of corresponding operators of 
the order A^^ j.} <^z,eS, which change the states of the 
resonators but not that of the qubit. 



IV. NUMERICAL RESULTS 

We now investigate the validity of the effective Hamil- 
tonian (|28p with respect to the resonator-resonator cou- 
pling constant, Eq. ()29p and the effective quantum master 
equation ([M)) for the resonators. Therefore, we compare 
the analytical results derived in the previous sections 
to numerically exact results obtained with the Bloch- 
Rcdficld quantum master equation, Eq. (jH]) using the full 
Hamiltonian ([5]). For further convenience, we assume 
uniform resonator decay rates, k = ka = kb- In our 
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FIG. 3: (Color online) (a) Rabi oscillations between 
the resonators at the effective interaction strength 
ffGSw/27r = 0.00315 GHz (numerical value, analytically 
we find gGSw/27r — 0.00328 GHz): Occupation num- 
bers of resonator A (solid red lines) and B (dashed 
blue lines). The initial state is |(7)qb|l)A|0)s. (b) 
Time evolution of the qubit state in terms of (ctz) 
(solid black lines). The deviation of the ideal case 
(a,) — —1 always lies below 0.5%. Numerical parame- 
ters: n/2Tv = 3.5 GHz, Sq/2tt = 4 GHz, e/27r = -6.37 GHz, 
g/2n = 0.24 GHz, G/27r = 0.0022 GHz, T = 20 mK. 
Decay rates: Ja(^)/2tv = Jb{^)/2tt = 0.00035 GHz, 
J^{n)/2Tv = J,{Q.)/2-K = 0.035 GHz. 

numerical simulations we use conservative estimates for 
the qubit decay rates. This is to stress the effect of the 
qubit dissipation channel on the resonators. 

A. Rabi oscillations 

The observation of Rabi oscillations between the two 
resonators is a first feasible application to probe the two- 
resonator setup. A system prepared in the product state 
l5)qb|l)A|0)B is subject to a periodic exchange of the ex- 
citation between the resonators as long as their coupling 
is finite, ggw 7^ 0. The corresponding oscillation time is 
Tkabi = T^/gsw- The initial excitation could be provided 
to one of the resonators by means of an ancilla qubit. 
For this purpose, suitable protocols have recently been 
proposedj ^^i^^" — 

Figure [3lja) depicts the according behavior of the res- 
onator populations as a function of time. We find the nu- 
merically observed oscillation period to be in good agree- 
ment to iRabi- Note that we have already incorporated 
the effects of the dissipative environments modeled by 
the Bloch-Redficld master equation. Eq. (jH]) (see discus- 
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sions in the following sections). The time evolution of 
the qubit population (tr^) is plotted in Fig. [HJb). From 
this we can verify that the qubit remains in its ground 
state after weak initial transients. These findings sub- 
stantiate the validity of the effective Hamiltonian (p8| in 
the dispersive regime. 



B. Decay rates 

In the following wc are interested in understanding 
quantitatively the influence of the reservoirs on the two- 
resonator setup. For this purpose we first make an an- 
alytical estimation based on the effective quantum mas- 
ter equation ()34|) . which arc compared then to numerical 
results obtained with Eq. We investigate the time 
evolution of particular observables, the associated oper- 
ators of which are constants of motion with respect to 
the dynamics of the closed system. Thus, any dynam- 
ics is produced by the dissipators of the quantum master 
equations (|9]) or (|34p . respectively. 

At this point we recall that the effective Hamiltonian 
simply describes a set of two coupled harmonic oscillators 
as long as the qubit remains in its ground state. They 
are each coupled to independent noise channels, as well 
as to a joint channel of qubit-induced correlated noise via 
their "center of mass coordinate". The latter is defined 
as the bosonic operator ^4+ = {a + b)/ \pl. In addition 
we define the "relative coordinate" A- = {a — b)/ \pl. In 
terms of these normal modes the oscillators are not cou- 
pled. The associated number operators are constants of 
motion, ["Heff, ^±^±] = 0. For further progress we com- 
pute the time evolution of the averages {J^j^Aj^ using the 
effective quantum master equation for the two-resonator 
system ([M)) . Here we note that the evolution of any op- 
erator O without explicit time dependence is described 
by the adjoint of the quantum master equation ([M]). 



TV CO 



(36) 



with 



+ k{D^ [a]0) + k{D^ [b]0) + K^b{D^ [a + b]0) , 



dt{0) =tv{dtp( 



The adjoint Lindblad super-operators D'^ act on the op- 
erator O, according to 

D'^lXp = X'^OX -]^[X'^X,0]+, (37) 

cf. Eq. (1331). Evaluating this relation for the normal 
mode number operators O = ^±^± yields 



{A\A+){t) = (4A+)t=o e-^-+^^--'>' (38) 



and 
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FIG. 4: (Color online) Decay of the quantum switch out 
of the initial state |(7)qb|l)A|0)B. Numerical data obtained 
with the Bloch-Redfield QME (O and = 3 states in each 
resonator: Decay of the expectation values of the number 
operators corresponding to the "center of mass coordinate", 
{A^j^Aj^) (red solid lines) and "relative coordinate", (^1j4_) 
(blue sohd lines), compared to analytical estimates, Eqs. (|38p 
and (|39p (black dashed and dash-dotted lines, respectively). 
Decay rates are according to Eq. (|34[) . Parameters: see Fig. [3] 



The normal modes are thus expected to decay exponen- 
tially. Remarkably, these decays should occur at different 
rates. The qubit-induced noise channel only couples to 
the center of mass, which suffers enhanced decay. This 
becomes manifest in the contribution 2Kqb to the expo- 
nent in Eq. with the rate Kqb from Eq. ([55)) . The 
relative coordinate is not affected by the qubit noise chan- 
nel, however, and simply decays with the resonator decay 
rate k, see Eq. ([39)) . Formally, this is because of 

D'^[a + b]A'^_A_ =0. 

In order to test these analytical estimations based on 
Eq. ([34)) . we consider a decay scenario with the resonators 
initially prepared in the Fock states \1)a (resonator A) 
and |0)b (resonator B). The qubit is prepared in its 
ground state |(?)qb- We calculate numerically the time 
evolution of the number operators related to the "cen- 
ter of mass", {Al^j^Aj^), and the "relative coordinate", 
(A^lA_), and compare the decay characteristics to the 
ones suggested by Eqs. ([55)) and ([5^ . respectively. The 
results are depicted in Fig. for a particular set of pa- 
rameters. We find an excellent agreement of theory and 
numerical data. While {A!^_AJ) decays at a rate k, the 

decay of {Al_^A^) is enhanced by the qubit noise channel, 
resulting in a decay rate The latter finding is 

affirmed in Fig. jH^b), where we compare the analytical 
expression for the decay rate k + 2Kqb to correspond- 
ing numerical values that arc extracted from simulations 
of a decay scenario according to Eq. ([55)) . The qubit- 
induced decay rate Kqb is given by Eq. ([55)) . Here, we 
have chosen an explicit dependence on the qubit bias en- 
ergy e, which is adjustable in realistic experimental sce- 
narios, while all other parameters are usually fixed. As 
indicated by Fig. |5t^a), the dispersive description can be 
considered as valid for A{A,E,n} ^ 0.1. These findings 
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FIG. 6: (Color online) Decay of the quantum switch out of the 
initial coherent state \a) with a = 1. Simulations were run 
with A'^ = 6 states in each resonator. Numerical data: occupa- 
tions of the "center of mass coordinate", {A'^A+) (red solid 
lines, oscillating curve) and "relative coordinate", {A^_A-) 
(blue solid smooth line), compared to analytical estimates, 
Eqs. (|38p and (|39|) (black dashed and dash-dotted lines, re- 
spectively). Initial transient effects are not depicted in the 
figure. Parameters: see Fig. [S] 



FIG. 5: (Color online) Dependence of the parameter Aa (a) 
and the effective damping rate K + 2Kqh (b) (black solid lines), 
on the qubit energy bias e. The rate for qubit- enhanced 
decay Kqb is given by Eq. (|35|) . The numerically extracted 
damping rates (red crosses in (b)) are related to the decay 
of the expectation value of the number operator correspond- 
ing to the "center of mass", {A\^A+), out of the initial state 
l5)qb|l)A|0)B at different e. The dotted lines mark the limit 
of validity of the dispersive theory (see text). Parameters are 
chosen as in Fig. O 



number in the resonator. This scenario is mainly mo- 
tivated by experiment, where a coherent state in a res- 
onator can easily be prepared via a resonant drive. To 
investigate the decay behavior of the "center of mass" 
and "relative coordinate" for this scenario, we choose 
the initial state |5)qb|a = ^)a\P — Q)b- As depicted 
in Fig. [HI the predictions of the effective quantum master 
equation are again found to be in good agreement with 
our numerical simulations, apart from transient effects. 



suggest that the effective quantum master equation for 
the two coupled resonators describes the dissipative 
system behavior adequately in the dispersive limit. 

Furthermore, the stationary value of {A^_^_A_^_) is found 
to be different from zero. This stems from a static en- 
ergetic shift of the oscillator potential minima due to 
a small "effective force". The latter arises from the 
resonator-qubit coupling component r2Ao(a + a^)o'z in 
the full system Hamiltonian ([5]), which has been elimi- 
nated in the effective Hamiltonian (jC17|) by the transfor- 
mation The dependence between the original and 
transformed oscillator creators and annihilators can be 
expressed as a*-^-* — > a^^-' — Aji. First of all, this shift 
explains strong oscillations with frequency D, during the 
time evolution of the resonator states. In Fig. SI their ef- 
fect is visible, however they are not resolved. The equilib- 
rium value of the resonator "center of mass" population 
is shifted according to 

(AU^) ^ - 4X1 ■ (40) 

Following the same reasoning, the equilibrium value of 
{A^_A_) is zero. 

As a second example, we consider the case of each res- 
onator prepared in a coherent or Glauber state, la) = 
g-\a\ J2'^^ij-^\n) with |ap being the average photon 



C. Decay of entanglement 

The generation of entangled two-resonator states is a 
key application of the quantum switch. For this purpose, 
we recall the switching property of the two-resonator 
setup mentioned in Sec. IIII Bl that is, the possibility to 
switch on and off the effective coupling between the res- 
onators by balancing the coupling coefficient gsw given 
in Eq. While a similar approach to create entan- 

glement between two resonators based on Landau-Zener 
sweeps^ has been previously discussed in Ref. HH, we fo- 
cus on the following, suitable protocol: A finite interac- 
tion strength is initialized by tuning the qubit energy 
flux appropriately. After preparing the initial product 
state |(7)qb|l)A|0)s the two-resonator state |^/')(.av evolves 
according to 

|V^W)cav =cos(gg^i)|l)A|0)B + isin(gg^t)|0)^|l)B . 

(41) 

After a time t = Ton has elapsed, gsw is balanced back 
to zero. During the whole procedure the qubit remains 
in its ground state |5)qb and does not get entangled with 
the resonators. In particular. Ton = T^/^gsw results in 
the entangled two-resonator state \^\.)cav = (|1)a|0)_b + 
*|0)a|1)b)/\/2, whereas Ton = 37r/4(7g^ yields the state 
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FIG. 7: (Color online) Decay of concurrence C+ and C_ 
for the entangled initial states and |^-), (solid red 

lines, thick curve, and solid blue lines, thin curve), respec- 
tively. The switch setting off-condition gsw = is fulfilled by 
e/2TT = -8.916 GHz/27r. Here we use G = 0.00478 GHz, all 
other parameters are as in Fig. [S] The exponential decay cor- 
responds to Eqs. I|44p and (|45|) (black dashed and dash-dotted 
lines, respectively). 

I^Dcav = {\1)a\0)b - i\0)A\l)B)/V2. A photon trans- 
fer from one resonator to the other is accomplished with 

Ton = 7r/2g§I^V- 

In the above discussions we have disregarded decoher- 
ence for reasons of clarification. In realistic scenarios, 
however, dissipation and dephasing are present even in 
the case of short times Ton, which prevents the creation 
of perfectly entangled states according to the above de- 
scribed protocol. Beyond that, two- resonator entangle- 
ment once created, will decay with time, according to the 
effective two-resonator QME In this context, excess 
phase noise in the resonators may cause further adverse 
effects, which are not considered here (cf. Sec. Ill Ap . 
Thus, it is important to reveal the decay characteristics 
of particular entangled states that could be created in 
the two-resonator setup up to a good degree via specific 
switch setting protocols. For this purpose, we first fo- 
cus on the decay characteristics of the initial entangled 
two-resonator Bell states 

IV'±)cav = ^(|l)A|0)B±|0)^|l)B). (42) 

To quantify the entanglement we first assume that all dy- 
namics is restricted to the subspace {|00), |01), |10), |11)}. 
Thus, we face the dynamics of entanglement between two 
two-level systems. In this case, the concurrence C repre- 
sents an adequate measure of entanglement, given by^ 
C = max{^i — ^2 — £,3 — ^4,0}. The parameters de- 
note the ordered square roots of the eigenvalues of the 
matrix Pca.v{o-ya^)p*^^{a^a^) with pcav being the re- 
duced density matrix of the two-resonator state, and A 
and B labeling the respective resonator Hilbert spaces. 
This representation of the concurrence is quite general 
and suitable for numerical investigation. However, for 
the initial states |'0±)cav and linear superpositions hereof, 
one can obtain analytical expressions for the decay char- 
acteristics of the concurrence with the help of the ef- 



fective quantum master equation (p4)) . Since the only 
nonzero elements of the associated density matrices dur- 
ing the whole time evolution are poo, Pii, P22 and pi2, P21 
in the basis {|00), |01), |10), |11)} the concurrence is sim- 
ply given by 

Cit)^2\pMt)\. (43) 

It turns out that the decay characteristics of the den- 
sity matrix element pi2 depend on the initial two- 
resonator state. In particular, the time evolution of 
the concurrences C± (t) for the initial density operators 
(|^±)(V'±l)cav is found as 

C+{t) = e-'''+2'''>'')* (44) 
C_(t) = e-«*. (45) 

The reason for this particular behavior is that the state 
|V'-)cav lies in a decoherence-free subspace with respect 
to the dissipator D[a + b]. Thus, it is a robust state in 
the sense that it does not couple to the qubit-induced 
correlated noise sourcc i^^'^^ This statement is equivalent 
to the relation D[a + 6](|V'-)(^/'-|)cav = 0. On the con- 
trary, the initial state |'0+)cav is fragile in this respect, 
since T'[a-h&](|^+)(?A+|)cav 7^ 0. 

In Fig. [71 we compare the numerically calculated time 
evolution of the concurrence to the analytical results of 
Eqs. p4)) and (|45|) . finding good agreement. While the 
decay of C+(t) is enhanced due to the qubit dissipation 
channel, the time evolution of C-(t) is determined by 
resonator dissipation only (cf. Fig. ^ , in analogy to the 
findings of Eqs. ([55)) and ([M)) . We note that a corre- 
sponding behavior has been reported for correlated states 
of a chain of coupled qubits interacting with a common 
bathi^ The numerical result is found to be shifted with 
respect to the analytical curves, since other elements of 
the density operator, e. g. P33 become populated as well 
during the time evolution of the system state. This non- 
RWA feature stems from the full numerical treatment 
using the system Hamiltonian Eq. ([T]). 

These findings can now be employed to characterize 
the decay of entanglement for the initial states |'0±)cav 
For this purpose we express them as linear superpositions 
of the Bell states |-0±)cav [Eq. P^ ]. 

\i''±U.^l({l + ^M+) + (l±^M-)) . (46) 

Z \ / cav 

Consistently, we find that the analogously defined con- 
currences C^- can be expressed as a sum of the concur- 
rences of the initial Bell states [Eqs. ()44p and ([45]) ] as 

exit) = Ci(t) = i ^^-{^+2n^i,)t _^ g-«t^ (47) 

This has some interesting consequences. For short times, 
the decay out of both initial states [■0j_)cav is merely gov- 
erned by qubit-enhanced decay at a rate n + 2Kqh- In 
the limit of long times, however, one finds pure resonator 
decay at a rate k. We have confirmed this numerically in 
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FIG. 8: (Color online) Decay of concurrence C+ for the en- 
tangled initial state l^*!!.). The switch setting off-condition 
gsw = is fulfilled by e/27r — —8.916 GHz. Here we use 
G = 0.00478 GHz, all other parameters as in Fig.|3l The two- 
mode exponential (solid black line) corresponds to Eq. (|47|) . 
The additional exponentials are given by Eqs. (|44|l and (|45|) 
and highlight the different decay regimes. 

Fig. |5]by means of the concurrence C^(i) related to the 
initial state |'0+)cav 

In summary, we point out that it is possible to under- 
stand the time evolution characteristics of the entangle- 
ment in the system on the basis of the effective master 
equation (p4)) . We emphasize that the qubit- induced dis- 
sipation channel plays a crucial, selective role for different 
classes of initially entangled states. 
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FIG. 9: (Color online) Time evolution of the auto- and cross- 
correlations of the initial state |5)qb|Q ~ 1)a\I3 — 0)b' Decay 
of {XI} + {X%) + 2(XaXb) (upper part, red lines) and {Xl) + 
(Xg) — 2{XaXb) (lower part, blue lines). The exponential 
decay is compared to the analytical estimates given by (|38p 
and (|39|l (black solid lines). Initial transient effects are not 
depicted in the figure. For parameters and initial conditions, 
see Fig. [3] and Fig. [S] respectively. 



V. EXTRACTING DAMPING CONSTANTS BY 
(CROSS-) CORRELATIONS 

In the two resonator setup, it is possible to measure 
correlations and cross-correlations in terms of the expec- 
tation values ((a-ha^)^), {{b + b'^y) and {{a + a^){b + b'^)) 
with present techniques <^ In the following we propose a 
method how to extract the relaxation rates k and Kqb out 
of correlation measurements of such type. We define the 
oscillator "positions" Xa ~ a + and Xb = b + b^ . In 
analogy to Sec. IIVBI we find analytically that the quan- 
tity (Xl) = {XD + {X%) - 2{XaXb) decays at the rate 
K. For {XD = {XD + {XD + 2{XaXb) we find a decay 
with K -|- 2Kqb. In Fig. [5J we numerically substantiate 
these findings for the example of a coherent initial state 
l3)qb|a = = 0)s, discarding again transient effects. 

Thus, we point out the possibility to extract the res- 
onator decay rate k, by measuring the decay of {X"^). 
This allows in turn for deducing the rate Kqb related to 
qubit-enhanced decay by measuring the quantity {X'j^). 
From the latter measurement, it is further possible to 
gain information about the relaxation and dephasing 
rates of the qubit, 7(wqb) [Eq. p^ ] and 7^(0; — >■ 0), 
[Eq. p^ ]. provided that the system frequencies and 
resonator-qubit interaction strengths are known. More 
details about a possible experimental realization of such 
correlation measurements can be found in Refs. [gsI 
and[66l 



VI. CONCLUSIONS 

We have investigated a two-resonator circuit QED 
setup. In the dispersive regime, i. e. if the resonator- 
qubit detuning is much larger than their mutual coupling, 
it is possible to extract the relevant system dynamics 
by applying the unitary transformation (|23|) to the sys- 
tem Hamiltonian ([5]). The resulting effective Hamilto- 
nian (|28|) reveals that the qubit gives rise to a switchable 
coupling between the resonators via virtual excitations. 
This dynamical coupling adds to the direct resonator- 
resonator coupling. Balancing both contributions, the 
resonator-resonator interaction can be set to zero. Such 
a qubit-mediated interaction provides a physical realiza- 
tion for a quantum switch between the resonators. 

As a principal point, we have focused on the dissipative 
system properties that stem from the interaction with 
different environments. For weak system-bath coupling, 
it is possible to cast the time evolution of the reduced 
system state into a quantum master equation of Bloch- 
Redfield form, Eq. It is usually derived starting from 
the total microscopic system-bath Hamiltonian, Eq. 
Its character being quite general, it only offers limited an- 
alytic insight. To study the dissipative dynamics in the 
dispersive regime, it is preferable to obtain a more use- 
ful, effective analytical description of the dissipative sys- 
tem dynamics. To this end, we have applied the unitary 
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FIG. 10: (Color online) Sketch of the relevant decay mecha- 
nisms affecting the quantum switcli. The qubit induces ex- 
tra decoherence in liigher dispersive order to the effective 
resonator-resonator dynamics. The arrows mark the deco- 
herence channels, labeled by the corresponding decay rates 
H{A,B,qh}- Decay processes of fourth order in the parameters 
A{A,E,Q} are not displayed individually. 



Appendix A: Energy relaxation and pure dephasing 
rates of the qubit 

In this section, we derive a quantum master equation 
for the qubit alone, which allows us to identify the energy 
relaxation and pure dephasing rates of the qubit. Con- 
sidering only a qubit coupled to individual environments 
along the a'^ and cr^-axes in the laboratory frame, the 
qubit-bath Hamiltonian reads 



n 



transformation Eq. ()23|) to the total system-bath Hamil- 
tonian ([8]) and analogously obtained the transformed, ef- 
fective system-bath coupling operators. Applying stan- 
dard methods, we have derived a rather complex effective 
quantum master equation for the system state in the ro- 
tated frame. It can be simplified, however, assuming low 
temperatures and recalling that the qubit state does not 
change, and only the dynamics of the two-resonator sys- 
tem are of relevance. By tracing out the qubit degrees of 
freedom one arrives at the Lindblad-typc quantum mas- 
ter equation for the reduced two-resonator state. 

As a main result of this paper, we have found that 
the qubit induces decoherence on the resonator-resonator 
system via an additional noise channel that acts on the 
"center of mass coordinate" of the resonators. This effect 
stems from qubit energy relaxation and is of second order 
in the small dispersive parameters A{a.s}i whereas pure 
qubit dephasing only enters in fourth order. This result 
anticipates that the operation of the quantum switch is 
robust against low-frequency noise in the two level sys- 
tem. For reasons of clearness, our findings are again illus- 
trated in Fig.[Tni We have verified our analytical findings 
by numerical calculations, where we have taken into ac- 
count the full dissipative dynamics according to Eq. 
As detailed in Sec. IIVI by means of several examples, we 
have found an excellent agreement of the presented dis- 
persive theory and the numerical results. Here, in par- 
ticular, we have investigated the validity of the obtained 
resonator relaxation rates. With regard to generating 
entangled states, which is a key application of the quan- 
tum switch, we have examined the decay mechanisms for 
different entangled initial resonator-resonator states. 



E EhM4a 



fJ.=X,Z J 



(Al) 



where H'^^ = {he/2)a'^ + {h5Q/2)al^. Applying Eq. Q, 
we obtain the diagonal qubit Hamiltonian = 
{hwc^h/'2)az and also the qubit-bath Hamiltonian, 



■Htot,qb = "Hqb + E E ^3 

+ (sin Oaz + cos 9 Ox) E (^j. 
-I- (cos a^ — 



' 3,3: 



(A2) 



Starting from the density matrix Ptot.nh associated with 
Htot,qb and following the lines of Ref. [SJ, the Lindblad 
quantum master equation for the reduced qubit density 
operator gqb = Trbath[£'tot,qb] can be derived. To this 
end, the spectral decompositions o'-(w) = S{ujq\, — w)(T_, 
(T+(a;) = (5(a;qb 4- w)(T+, and crz(j-L>) — (5(a;)o'2 of the qubit- 
bath coupling operators are required. Omitting the ex- 
plicit time dependence of gqb for simplicity, we find 



eqb - - ^ 



- + h + - 



+ 7(^^qb) 
+ 7(-kqb|) 

+ 70(w->O) cr z Qqhcr z - gqh 



• 9qh\ 



o-^QqhCr' 



' fi'qb] + 



(A3) 
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Here, (T_ = \g){e\ and a+ = \e){g\ are the fermionic qubit 
annihilation and creation operators. The energy level 
transition and the pure dephasing rates are given by 

-i{uj)=^T^{uj)cos^e + T^{uj)siT?d (A4) 
7^(0;) = r^(a;) sin^ 9 + T ^{uj) cos^ 9 , (A5) 

respectively, and depend on the bath correlation func- 
tions 



J^(|w|)n^(|a;|) w<0. 



(A6) 



where n^{uj) = ^^^^gRw/feeT,, _ 2) is the Bosc distribu- 
tion function of bath with label /i S {a;, z}. Because 
the quantum switch operates in the limit of low tem- 
peratures, k^Tfi <ti fixjqb, the Bose-factor van- 
ishes for frequencies of the order of Wqb- However, for 
a; — )• 0, n^{ijj) tends to diverge. This can be relevant 
in the experimentally important case of l/f-noise?^"— 
which would require a treatment beyond the framework 
of a Markovian master equation, exceeding the scope of 
this work. Instead, we avoid the divergence problem by 
choosing Ohmic spectral densities [Eq. pT|) ] for low fre- 
quencies w ^ Wqb. In many cases, this assumption is 
reasonable.—"— Provided that both baths have the same 



temperature T = Ti 



{x,z} 



we obtain 



'-f{uj > 0) = Jx{io) cos^ 9 + Jz{uj) sir? I 
7(w < 0) « 



70 



(oz cos 9 + sin 



(A7) 
(A8) 

(A9) 



These rates are functions of the mixing angle 9. Eq. (jA7[) 
constitutes the main result of this section. For u> = Wqb, 
it establishes the connection between the energy relax- 
ation rate in the energy eigenframe and the dissipative 
baths defined in the laboratory frame. The Markovian 
description of Eq. (|A3p remains justified as long as the 
spectral densities J{x.z}{i-^) are smooth functions in Wqb. 
This allows one to apply a linear approximation in an 



infinitely small interval around ojqb, which yields an ef- 
fectively Ohmic description. We finally note that in the 
special case of noise coupling purely via the laboratory 
z-axis (e.g., a flux qubit exposed to flux noise) the re- 
sult 7(a;qb) = Jz i'^qh) sin^ 9 is in agreement with findings 
from other works.— 1^>^ For the sake of completeness 
we mention that the master equation (jA3|) together with 
the rates (jA7p . (jASp and (|A9p reproduce the well-known 
results concerning relaxation and pure dephasing times, 
(Ti)-i = 7(0; > 0) and Tr^ = -f^{uj 0) respectively. 



Appendix B: Effective bath coupling operators 

To obtain the quantum master equation for the re- 
duced system state, it is necessary to transform the total 
system-bath Hamiltonian to the dispersive picture via 
the transformation Htot,off = U^Htod^, with h{ given 
in Eq. The effective system Hamiltonian having 

been derived in Sec. IIII Bl we need yet to find the trans- 
forms (j23|l of the system-bath coupling operators Qfj, to 
the dispersive frame. Up to second order in A{a,e,o}i 
they read as 

fi e {A,B,x,z} . 

(Bl) 



^(^{A,E,f2}) 



Each of the effective bath coupling operators Q^^cS is 
represented by of a sum of operators. 



For the resonator-bath coupling operators Qa = a + 
and Qb = b -\-b\ wc obtain the dispersive transforms as 



Qyi.eff = (a + a' )cff =a + aJ + (Aa - X^)(Tx - 2\n(yz 

+ ^(Aa - As')(a + h + a^ +b^) + Ao(Aa + AE)a,(a + 6 + + 6^) , (B2) 
Qb,cS = (& + fc^)cfr =5 + 6^ + (Aa - As)o-^ - 'iXncTz 

+ ^(Aa - As')(a + 5 + at + fet) + ^^(Aa + X^)a^{a + 6 + + 5^) . (B3) 

I 



The dispersive transforms of the qubit-bath coupling op- 
erators Qx,eS = <yx,cS and Qz,cS = Oz,cS are obtained as 
a combination of those of the qubit operators cr^ and cr^ 



in the laboratory basis, according to Eq. The latter 
assume the form 
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(^'x,cs + (Aa + As)cr^(a + 5 + + 6^) - 2Af2 a"^ {a + h) + {o) + 6^) - (T+(a^ + fe'f) - cr"(a + &) 

- ((Aa + Ae)2 - 4A^)(T, [(at + 5t)(a + 6) + l] + - (Aa + Ae)Aa) [(7^(0^ + 6^)' + (T+(a + 5)^] (B4) 
+ An(AA - As)(T,(a + 6 + at + 5t) + (2A2 - (Aa + As)As) [c7-(a + 5)^ + a^{a^ + &t)2 



eff =0-2 - 2Aa cr+(a + 6) + (t (at + fot) 



- 2Aj 



CT+(at + &t) + + 



(aI + A|) [(at + 6t)(a + 6) + l] + 4Ao(Aa - As)(t. [(at + 5t)(a + 6) + 1 



4AoAa 



(T+(a + 5)2 + cr-(at + &t)2 



+ 4Aj2As 



CT-(a + &)2 + (T+(at + 6t)5 



2AaAso-2 



(a + &)2 + (at + 6t)2 
(B5)' 



Appendix C: Effective quantum master equation in 
the dispersive limit 

Starting from the Bloch-Redficld quantum master 
equation we move to an interaction picture with re- 
spect to the system and the individual reservoirs. Here, 
the coupHng operators have to be replaced by their 
dispersive transforms Q^^off found in App. |B] Now, we 
introduce the spectral decompositions 



(Cl) 



The Qj^ are the summands of the effective coupling oper- 
ators as detailed in Eqs. (|B2[) - (jB5p . The spectral compo- 
nents Qj^ (cj) are obtained by expanding the Qj^ in terms 
of the eigenstatcs of the effective Hamiltonian which 
we cast in entirely diagonal form for this reason, 

U,^ = hh^ [a\A^ + i) (^t_ a_ + 1) + ^ , (C2) 

Here we have defined 

^+ = r2 + G + 2(AAA + AEl])cr2 (C3) 

0_ = O - G (C4) 
? = a;qb + 2(AAA + AsS) , (C5) 

and introduced via a linear transformation the normal 
modes 



A 



1 



1 



(C6) 



The eigenstates of the effective Hamiltonian (|C2p can be 
considered as re-defined Fock states. 



-I- h£L_ ( m + i ) |m) 



(C7) 



Here, {n,m} = {0,1,2,...} denote the oscillator exci- 
tations (or resonator photon numbers), and f2yi(0 can 
assume the values 



fi+(0 = f7 + G 2(AaA AsE)(-l 



(C8) 



with ^ = 1 or ^ = for the qubit being in the excited or 
ground state, respectively. With this at hand we find the 
spectral decompositions of the effective coupling opera- 
tors via the ansatz 



nml n'm'l' 

X \nml){nml\Qj \n'm!V) {n'm'l' 



(C9) 



where A"JJJ ' denotes the energy difference between the 
states \n'm'l') and \nml). For illustration we list the ex- 
plicit expressions for the spectral decompositions of some 
components. 



cr"(a;) 



A+|0)(0|5(c.-f}+(0)) 
+ A\l){l\5{Lo-n+{l)) 
A^5{uj-n^) 

\n){n\5[Lj- [wqb + (2n + l) 



El 



X (AiA 



■A|S)] 



(CIO) 
(Cll) 

(C12) 

(C13) 

t 



The decompositions of operator products such as A 
etc. are obtained analogously via the relation (w) = 
Qj^X^\io\). With this and Eq. (jCl|) . we recast the Bloch- 
Rcdfield quantum master equation ([9]) into the form 



X {^Q,^{u:)Q{t)Q\y)~Q\y)Q,^{uS)Q{t)) + h.c. , 

(C14) 
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where Ti^i{lu) is the one-sided Fourier transform 



/>oo 

r^(w) = / dre^'^^K^iT) 
Jo 



(C15) 



of the bath correlation function Ky^{T) given in Eq. ([10)). 
One usually neglects the Cauchy principal value of the 
integral and can then rewrite (jClSP as 



J^(Lj)(n^(w) + l) w>0 
Jp(|w|)n^(w) w<0, 



(C16) 



with the spectral density J^{oj) and the Bose distribution 
function n^{ijj) = \j{e^^l^'B.'^i^ — 1), depending on the 
temperatures T^. 

Inserting the explicit expressions for the spectral 
decompositions into the quantum master equation 
(|C14p . we find two different classes of oscillating 
terms. The first oscillate at high frequencies such as 
giifiA/B*^ gijt^qbt^ giiAt^ g±iEt^ i.e. vary on timescales 
of the intrinsic system evolution, whereas the second 
oscillate slowly at frequencies A^A + A|;S and multi- 
ples. This difference enables one to perform a semi- 
secular approximation similar to the approach in Ref. Issl . 
Here, we assume that all rapidly oscillating terms of the 
first class can be averaged to zero. This is justified 
since the timescales of intrinsic system evolution given 



by (f^A/s)""'^ etc. are typically much smaller than the 
relaxation timescales, on which the system state varies 
notably. This, however, is not the case for terms of the 
second class, which we keep consistently. We empha- 
size that our approach goes beyond the standard way of 
obtaining a Lindblad quantum master equation. The lat- 
ter would imply a full secular approximation, neglecting 
all oscillating contributions and only keeping terms with 
uj^J in Eq. ((CT4)) . 

Furthermore, we may simplify the bath correlation 
functions, 

r(wqb + n{XlA + AlS)) « FKb) 

for small occupation numbers n, and assume an over- 
all temperature T = Ti^a.b,x,z}- In the low-temperature 
regime, T <C (fi./fcB) minjwqb. A, 57}, it is appropriate 
to neglect all contributions to Eq. (|C14[) with negative 



frequencies because of F^(a; < 0) 



0, 



ergy is absorbed from the baths. This automatically 
yields F^(w) = J^{uj). In the low- frequency region of 
the qubit baths, we assume Ohmic spectral behavior, 
J{x,z}i^) = OL{x,z}'^- As detailed in App. 1X1 this implies 
^{x,z}{^ ^ 0) = a{^^^}kBT/h. 

We eventually obtain the effective quantum master 
equation for the reduced system state 



- [Hoff, p] + JA{n)D[a]p + JB{n)D[b]p + (Aa + \^f{Jx{^) cos' 9 + J,{n) sin' 9)D[cj,{a + b)\p 
+ J, {ujqh)D \a- ( cos - ( cos 9{-4Xl + (Aa + As)') + 4 sin 9Xn (Aa - As )) ((a^ + 6^ ) (a + 6) + l) ) 
+ J,(Ljqb)i^[(T-( - sine + (sin0(-4A^ + (Aa + As)') - 4cos0Ao(Aa - As)) ((a^ + 5t)(a + &) + l) ) 
-t- 4(J:^(A)(An cos6' + Aa sin 6*)' + J^(A)(Aa cos6I - Ao sin 6*)') Z) [(T-(a^ + fc^)j P 
4-4(j^(I])(Ancos6'-Assin6i)' + J2(I])(AsCos6' + Ansin6l)')i:) cr"(a + &) p 
+ {kBT/h){a^ sin' 9 + az cos' 9)D a, (l - 2(A^ + A|) ((a^ + &^)(a + &) + l)) 

I 



(C17) 



where we have omitted the time dependence of the den- 
sity operator p and used the notation D[X]p = XgX^ — 
i [X^'X, g] ^ with the anti-commutator [A, B]+ — AB + 
BA. 

We have neglected terms cx -^'a s f2}'^{'4--B} ('^) ^"^'^ 
higher orders of AjA.s.n} in (|C17p . This is justified, since 
the resonators typically employed in experiment possess 
a high quality factor and therefore feature small decay 
rates. On the other hand, it is necessary to keep terms 
proportional to A'^ ^ n}'^{a;,z}(^)' because typical qubit 
dephasing and relaxation rates exceed those of the res- 
onators by several orders of magnitude with conserva- 



tive estimates. In general, it is appropriate to neglect all 
terms of the order of A|^ ^ and higher. Due to the dis- 

sipator relation D[AX] = A'Z3[X] we thus may already 
discard any contributions to the system-bath coupling 
operators of order A'^ ^ in Eqs. (jB2p - (|B5|) . 



In Sec. IIII Al we have motivated the experimental ad- 
vantage of keeping the qubit in its ground state without 
own dynamics. This enables us to simplify (jC17p further 
by tracing out the qubit degrees of freedom. To this end, 
we take into account "T'\:{(Tz\g){g'\) = —1 and the partial 



17 



trace relation 

Tri(D[Ai ® A2]p) = Tri(AlAiD[A2]p) , (C18) 

with Ai and A2 acting each on a different Hilbert sub- 
space, and Tri(-) denoting the partial trace with re- 



spect to one subspace. Finally, we arrive at the effective 
quantum master equation for the tow-resonator state, 
Eq. ([34|) . The qubit decoherence rates and 7^ are 
identified following the discussion in App. [X] 
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